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I.  INTRODUCTION 

The  rapid  growth  of  offshore  oil  production  and  marine  oil  transport 
has  led  to  an  increasing  danger  of  oil  contamination  of  the  coastal  environment. 
Consequently,  there  is  a  growing  worldwide  concern  over  possible  environ¬ 
mental  damage  caused  by  accidental  oil  spills. 

The  present  study  has  been  aimed  at  developing  a  computer  model 
which  is  capable  of  predicting  oil  slick  transport  in  harbors  and  bays.  By 
means  of  numerical  approximation,  the  model  is  able  to  simulate  the  spread¬ 
ing  and  the  movement  of  oil  slicks  on  the  ocean  surface  and  to  predict  their 
destinations. 

The  movement  of  an  oil  slick  is  simulated  as  a  combination  of  the 
various  phenomena  which  affect  the  spreading  and  transport.  The  simulation 
of  the  spreading  process  is  governed  by  the  physical  properties  of  oil  and  its 
characteristics  at  the  air-oil  and  oil-water  interfaces.  The  spreading  motion 
is  then  superimposed  on  the  drift  motion  caused  by  winds  and  tidal  currents 
to  give  the  total  movement.  By  considering  the  slick  as  a  summation  of  many 
elementary  patches  in  its  numerical  scheme  and  applying  the  principle  of 
superposition  to  each  individual  patch,  the  model  is  capable  of  simulating  the 
slick's  shape  distortion  as  a  result  of  the  relative  shear  motion  produced  by 
the  non- uniformity  of  tidal  currents  and  wind  drifts  on  the  water  surface. 

Sample  computations  have  been  conducted  using  the  Long  Beach  Harbor- 
San  Pedro  Bay  area  as  a  test  case.  Several  field  tests  using  either  cardboard 
markers  or  soybean  oil  to  simulate  a  spill  were  also  conducted  in  this  area. 
Computer  predictions  using  input  obtained  from  field  information  show  good 
agreement  with  the  field  traces  and  the  validity  of  the  computer  model  for 
field  operation  is  substantially  confirmed. 

In  order  to  accommodate  the  computer  model  for  use  in  local  port 
offices,  a  dedicated  small  computer  system  and  an  alternate  time- sharing 
system  are  proposed.  In  either  case,  the  geographical  data  and  the 
oceanographic  environment  of  a  harbor  can  be  pre-fed  into  the  system. 
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Wind  information  may  be  fed  constantly  to  the  system  from  remote  sensors 
in  the  field.  In  the  event  of  an  accidental  spill  occurring  in  this  area,  one  may 
press  a  keyboard  to  supply  very  simple  data  such  as  the  time  of  the  spill, 
type  of  oil,  location  and  volume  of  the  spill,  and  from  this  input  the  model 
will  predict  the  oil  size,  shape,  and  movement  as  a  function  of  time.  The 
model  may  provide  a  moment-to- moment  update  of  critical  conditions.  The 
model  may  predict  the  arrival  time  of  the  oil  slick  at  various  portions  of 
the  shoreline  along  with  its  spreading  and  transport  well  ahead  of  the  true 
slick's  arrival.  This  advance  information  of  arrival  time,  location,  and 
size  would  be  extremely  useful  in  guiding  the  decisions  of  engineers  and 
operational  personnel  for  an  effective  response  in  control  and  cleanup, 
and  therefore,  the  model  can  be  a  very  effective  tool  to  minimize  the 
environmental  impact  on  shorelines  as  well  as  on  the  water  quality  in 
harbors . 


II. 


OIL  SPREADING  ON  CALM  WATER  SURFACE 


Oil  spills  generally  occur  in  two  forms.  They  can  be  in  the  form 
of  a  continuous  release  with  a  time-dependent  spill  rate  such  as  that  which 
occurred  in  the  Santa  Barbara  Channel  as  a  result  of  terminal  failure,  or 
they  can  be  in  the  form  of  a  nearly  sudden  release  such  as  which  may  occur 
as  a  result  of  a  marine  accident  of  oil  carriers.  The  distinction  between 
the  two  is  mainly  dominated  by  the  volume  of  the  spill  involved.  The  failure 
of  an  oil  drilling  terminal  naturally  involves  a  large,  volume  and  long  time 
and  should  be  characterized  as  a  continuous  spill.  On  the  other  hand,  a 
relatively  small  volume  spill  due  to  an  accidental  dumping,  which  although 
can  hardly  complete  instantaneously,  may  generally  be  regarded  as  a 
sudden  release. 

In  this  report,  only  a  sudden  release  of  oil  is  considered.  The  site 
of  spill  is  considered  inside  a  harbor  or  inside  the  breakwater  protection  of 
a  bay  where  the  effect  of  the  waves  are  generally  negligible  and  the  water 
may  be  regarded  as  calm. 


The  basic  mechanism  affecting  the  spreading  of  oil  slicks  over 
calm  water  includes  inertia,  gravity,  viscosity  and  surface  tension.  At 
the  initial  stage,  the  primary  driving  force  is  due  to  gravity  and  the 
rate  of  spreading  is  governed  by  a  balance  between  the  gravitational 
pressure  and  the  oil  inertia.  As  the  spreading  proceeds,  the  oil  slick 
becomes  thinner  and  the  viscous  effect  becomes  more  evident.  Quickly, 
the  oil  thickness  becomes  thin  enough  so  that  the  inertia  effect  becomes 
negligible  and  the  spreading  enters  into  its  second  phase  while  the  gravita¬ 
tional  spreading  force  is  primarily  balanced  by  the  water's  viscous 
retarding.  The  gravity  oriented  body  force  gradually  becomes  less 
important  as  compared  to  the  air-water-oil  interfacial  effect  when  the 
film  thickness  becomes  even  thinner,  typically  on  the  order  of  mili- 
meters.  Finally,  surface  tension  becomes  dominant  as  the  major  driving 
mechanism  and  responsible  for  the  final  phase  of  spreading.  The  spread¬ 
ing  rate  at  this  phase  is  therefore  determined  by  a  balance  between  the 
surface  tension  spreading  force  and  the  viscous  drag. 


By  balancing  the  various  combination  of  forces,  Fay  [)]  derived 
the  governing  equations  for  various  regimes.  The  actual  spreading 
processes  at  various  regimes  were  also  investigated  experimentally  and 
results  were  compared  with  analytical  predictions  (Moult  [2]).  While  the 
experiments  show  that  the  theory  predicts  correctly  in  terms  of  the 
functional  relationship  among  the  governing  parameters,  discrepancies 
exist  on  the  magnitude  of  the  proportionality  constants.  Waldman  et  al  [3], 
also  recommended  their  values  of  these  proportionality  constants  based 
upon  an  analytical  prediction  pc  formed  by  Fannelop  and  Waldman  [4], 
Their  values  in  some  cases  cuuer  as  much  as  50%  from  Fay's  recommen¬ 
dation.  In  the  following,  the  governing  equations  for  radial  spreading 
derived  by  Fay  are  summarized.  With  understanding  that  large  discrep¬ 
ancies  exist  between  theories  and  experiments,  the  proportionality  con¬ 
stants  were  chosen  mainly  based  upon  fitting  of  empirical  data. 
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spreading  coefficient  (net  surface  tension 
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For  a  small  volume  spill,  only  the  third  phase  is  of  significance 
as  the  first  two  phases  last  only  a  very  short  period  in  the  entire  spread¬ 
ing  history.  For  instance  of  a  10,000-gallon  spill,  the  inertial  phase 
lasts  only  20  minutes  and  the  viscous  spreading  process  should  complete 
in  40  minutes.  Consequently,  Equation  (3)  alone  is  sufficient  to  describe 
the  spreading  almost  entirely  except  in  the  very  first  hour  after  the  spill 
starts.  Good  agreement  of  Equation  (3)  with  the  observed  data  may 
show  the  appropriate  choice  of  the  proportionality  constant  (Figure  1).* 

Equation  (3)  also  shows  that  within  the  surface  tension  process, 
the  leading  edge  of  the  oil  spreading  is  independent  of  the  volume  spilled. 
There  is  no  question,  however,  that  the  final  area  of  the  film  does  depend 
on  the  volume  of  the  oil  spill.  Theoretical  determination  of  the  final  area 
would  involve  the  knowledge  of  various  physical  properties  of  oil  in  water, 
such  as  the  change  of  the  spreading  coefficient,  the  diffusivity,  and  the 
solubility.  For  practical  purpose,  however,  an  overall  estimate  proposed 
by  Fay  [5]  is  used;  the  relationship  is  given  as  follows: 

A  (m2)  =  105  [V(m3)]3/4  '  (4) 

The  expression  seems  much  over-simplified;  it  relates  with  only  one 
single  parameter,  the  spill  volume.  Nevertheless,  it  does  fit  the  field 
data  very  closely  as  shown  in  Figure  2. 


♦Figures  are  given  at  end  of  report 
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Nominal  values  for  tlic  physical  properties  of  water  and  crude 
oil  arc  used  in  the  entire  study;  they  arc  listed  as.  follows: 

3 

p  ,  =1  gm/cc  =  1.94  slug/ft 

rwater  b  b 

=  0.95  gm/cc  =  1.84  slug/ft^ 

V  =  0.  012  cm^/scc  =  1.290x10  ^  ft^/sec 

0  =  9.  5  dynes/cm  =  0.  65  x  10  ^  lb/ft 

Based  upon  these  nominal  values,  calculated  results  showing 
various  regime  durations  are  presented  in  Figure  3.  This  figure  clearly 
shows  that  in  small  volume  spills  the  surface  tension  phase  is  dominant 
during  the  entire  spreading  history. 

o  A 

Calculations  of  the  slick  growth  for  spill  volumes  of  10  ,  10  , 

10^,  and  10^*  gallons  are  presented  as  a  function  of  time  in  Figure  4. 
Since  the  slick  size  is  independent  of  volume  in  the  surface  tension 
phase  of  spreading,  all  patches  propagate  along  the  same  path  line 
until  they  reach  their  final  dimensions  and  cease  to  spread. 

The  thickness  of  the  oil  slick  varies  from  its  center  to  its  edge. 
The  exact  thickness  distribution  of  a  slick  can  only  be  obtained  through 
proper  modeling  of  the  spreading  dynamics  during  its  growth.  In  order 
to  give  an  order  of  magnitude  idea  about  the  slick  thickness  during  the 
various  phases  of  spreading,  however,  Figure  5  has  been  prepared  by 
assuming  uniform  thickness  distribution  for  all  cases  at  all  times.  The 
figure  shows  the  nominal  thickness  variation  as  a  function  of  time  for 
the  same  four  spills. 
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III. 


OIL  TRANSPORT  IN  HAP  BOR 


The  primary  factors  affecting  the  movement  of  an  oil  slick  are 
winds,  tidal  currents  and  waves.  In  the  foregoing  discussion  on  the 
spreading  process,  the  water  ‘surface  has  been  assumed  calm  without 
any  effect  due  to  water  circulation  or  perturbation.  In  describing  the 
transport  process,  the  assumption  is  simply  extended  that  the  drifting 
motion  of  an  oil  slick  caused  by  winds,  currents  and  waves  can  be  super¬ 
imposed  on  the  spreading  motion  of  the  slick  on  calm  water  as  described 
in  the  previous  section. 

The  direct  drift  caused  by  the  wind  shear  stress  over  t'.e  water  sur¬ 
face  is  generally  agreed  being  on  the  order  of  3%  of  the  over-water  wind 
speed.  There  is  disagreement  as  to  the  direction,  however.  For  in¬ 
stance  the  experiments  of  Teeson  and  Schenck  [6]  and  the  analysis  of 
Warner,  et  al  [7]  suggests  that  an  angular  displacement  on  the  surface 
drift  due  to  the  effect  of  earth  rotation  must  be  considered,  whereas  the 
analysis  of  the  Torrey  Canyon  Oil  Slick  Movement  [8]  shows  good  agree¬ 
ment  by  assuming  the  oil  drift  always  aligned  with  the  wind  direction. 

S 

It  is  anticipated  that  Coriolis  force  will  affect  the  course  of  trans¬ 
port  if  an  oil  slick  is  in  movement  under  a  steady  wind  over  a  large  distance 
on  an  open  sea.  While  inside  a  harbor  as  opposed  to  an  open  sea,  the  winds 
are  highly  local  and  unsteady  in  some  cases  and  the  Coriolis  effect  can 
hardly  be  developed  fully  in  a  short  time  scale;  nevertheless,  a  partial 
development  of  the  Coriolis  effect  may  be  significantly  responsible  for  oil 
transport. 

The  effects  of  surface  waves  are  generally  negligible  because  of 
their  oscillatory  nature  which  produces  little  net  force  affecting  the 
spreading  motion.  Certainly,  nonlinear  effects  of  wave  interaction  may 
produce  surface  drift  due  to  second  order  mass  transport.  On  the  other 
hand,  however,  there  is  a  retardation  effect  caused  by  waves  on  sur¬ 
face  drifts  because  of  the  presence  of  wind  shadow  on  the  lee  side 
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of  the  wave  crest.  Schwartzberg  [9]  found  that  the  enhancement  and  the 
retardation  in  surface  drift  caused  by  waves  are  approximately  of  the 
same  order  of  magnitude  and  cancel  each  other.  The  apparent  wind  drift 
can  factually  be  regarded  as  a  constant,  and  no  net  effect  caused  by  waves 
needs  to  be  considered. 

The  magnitude  and  direction  of  the  tidal  current  are  assumed,  in 
the  present  work,  either  calculable  or  measurable  with  required  accuracy. 
In  the  numerical  example  given  in  the  later  part  of  this  report,  the  tidal 
currents  were  calculated  as  a  function  of  time  and  space  with  the  tidal 
stage  as  the  input  data. 

In  summary,  the  transport  process  of  an  oil  slick  is  handled  as 

follows: 

1.  The  effects  of  wind  and  current  are  assumed  uncoupled 
with  the  spreading  dynamics  and  they  are  superimposable 
on  the  spreading  motion  of  an  oil  slick  on  calm  water. 

2.  The  effects  of  waves  on  surface  drift  arc  considi  red 
negligible.  '  - 

3.  The  tidal  currents  are  obtained  from  either  calculations 
•  or  measurements.  They  are  input  data  varying  as  a 

function  of  time  and  space. 

4.  The  wind  induced  surface  drift  is  assumed  to  have  a 
magnitude  of  3%  of  the  wind  speed.  Its  angle  of  deviation 
from  the  wind  direction  is  to  be  provided  as  a  part  of  the 
input  information;  otherwise  it  is  assumed  to  be  aligned 
in  the  same  directio  n  as  the  wind. 


IV. 


NUMERICAL  SIMULATION 


1.  Model  Set  Up 

The  two-dimensional  field  of  a  harbor  is  constructed  by 
a  set  of  perpendicular  grids  of  equal  grid  spacings.  At  each  grid 
point,  the  field  characteristic,  ie,  ,  land  or  water,  is  identified 
by  information  fed  from  input  data.  Each  grid  point  represents 
a  rectangular  area  of  As  x  As  ,  if  As  is  the  grid  spacing. 
Whenever  the  grid  point  represents  water  in  the  field,  the  local 
water  depth  should  also  he  specified.  In  the  two-dimensional 
field,  the  tidal  current  and  the  wind  driven  current  are  each 
represented  by  a  two-dimensional  velocity  vector  at  each  grid 
point.  The  variation  in  magnitude  and  direction  of  these  vectors 
from  one  grid  point  to  another  signifies  the  surface  current  dis¬ 
tributions.  Similarly,  if  there  were  oil  slicks  in  the  field,  the 
local  film  thickness  (averaging  over  a  As  square)  will  be 
specified  at  the  corresponding  grid  points. 

2.  Simulation  of  Oil  Spreading  and  Transport 

The  spreading  of  an  oil  slick  consists  of  three  distinct 
phases  as  discussed  in  Section  II.  The  relationships  governing 
the  spreading  phenomenon  in  different  phases  presented  in  Section  II 
provide  the  fundamental  information  of  circular  spreading  for  an 
oil  slick  on  a  calm  water  surface.  When  the  oil  slick  is  on  water 
of  uniform  current,  the  effect  of  the  water  movement  may  simply 
be  superimposed  to  the  spreading  motion.  When  the  film  is  subject 
to  non-uniform  disturbances  resulting  from  tide,  wind  and  others, 
the  slick  shape  would  be  distorted  and  appear  irregular.  To 
solve  the  exact  problem  concerning  spreading  and  transport 
in  an  arbitrary  velocity  field  requires  a  detail  analysis  of  forces 
on  elementary  slick  segments.  The  analysis  would  very  likely  re¬ 
sult  in  a  set  of  differential  equations  similar  to  the  convective 
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diffusion  typo  equations  with  a  tensor  of  air-oil -water  interface 
spreading,  equivalent  to  a  tensor  of  turbulent  mixing  relating 
mass  transfer  between  streamlines  in  the  pollutant  transport 
problem.  Assuming  that  this  spreading  coefficient  could  he 
specified  as  a  function  of  space  and  time,  there  should  be  no 
difficulty  in  general  to  write  a  computer  program  to  determine 
the  movement  and  spreading  of  an  oil  slick  in  a  given  environ¬ 
ment  as  well  as  the  slick  thickness  distribution  as  a  function  of 
space  and  time.  Nevertheless,  even  with  modern  computers  such  an 
attempt  would  be  too  expensive  to  be  practical  in  terms  of  computing 
time,  if  not  just  ambitious. 

Instead  of  solving  the  exact  problem,  the  present  model 
simulates  the  spreading  and  transport  phenomena  numerically  by 
means  of  superposition.  The  oil  film  is  fictitiously  divided  into 
patches  acc<  ding  to  the  already  established  two-dimensional 
grids;  each  patch  occupies  an  area  of  As  square  centered  at  the 
grid  point.  Th<  concept  is  to  follow  the  motion  of  the  center  of 
these  patches  and  to  superimpose  their  motion  together  with  a 
spreading.  The  process  is  carried  out  in  appropriate  time  intervals. 
At  the  end  of  each  time  step,  the  average  thickness  of  each  patch  is 
obtained  by  summing  up  the  overlapping  layers,  resulting  from 
spreading  and  transport  at  each  grid  point. 

While  the  detailed  thickness  distribution  is  not  exactly 
material  in  the  present  study,  the  development  of  the  slick  sire 
must  be  correctly  simulated.  It  has  been  discussed  in  the  fore¬ 
going  sections  that  the  slick  size  as  well  as  the  regime  of  spread¬ 
ing  depends  upon  the  total  volume  of  oil  and  the  total  time  elapsed 
from  the  beginning  of  a  spill.  Instead  of  following  the  total  volume 
of  the  oil,  the  present  scheme  focuses  its  attention  on  the  patches 
at  each  grid  point.  In  fact,  the  adopted  concept  considers  that 
each  patch  spreads  under  the  influence  of  its  own  volume.  At 
the  end  of  each  time  interval,  new  patches  of  new  volume  content 
are  to  be  established,  and  each  new  patch  spreads  accordingly 
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with  its  new  volume  in  the  next  time  interval.  Following  this 
concept,  at  the  beginning  of  each  time  step,  all  patches  are 
essentially  being  considered  as  new  spills  without  previous 
spreading  history  and  the  effective  time  of  spreading  for  each 
of  these  patches  therefore  counts  only  from  the  beginning  of  the 
time  interval  of  concern.  The  concept  is  similar  to  that  adopted 
by  Fischer  [10]  in  predicting  pollutant  transport  in  water.  As  a 
result  of  carrying  out  this  concept,  the  center  of  the  slick  is 

thicker  than  elsewhere  around  initially,  but  the  gradient  of  the 

\ 

film  surface  decreases  gradually  till  a  limiting  thickness  governed 
by  the  spreading  mechanism  is  reached. 


In  numerical  computations  dealing  with  the  diffusion  type 
equations,  the  time  step  and  grid  spacing  must  be  kept  in  appro¬ 
priate  limits,  depending  upon  the  magnitude  of  the  diffusion 
coefficient,  in  order  to  assure  the  solution  to  be  stable.  Similarly, 
in  the  present  numerical  scheme,  the  correct  simulation  of  the 
spreading  mechanism  requires  an  appropriate  choice  of  the  grid 
size  and  the  time  step.  The  selection  of  them  depends  mainly  on 
the  spreading  coefficient  and  the  spill  volume. 

It  is  understood  that  the  inertial  and  viscous  regimes  are 
generally  short  as  compared  to  the  entire  spreading  history. 

Since  the  model  is  intended  to  be  simple  enough  for  a  quick,  order 
of  magnitude  estimation,  in  most  cases  the  time  interval  and  the 
grid  space  can  be  chosen  in  such  a  way  that  only  the  last  phase 
(surface  tension)needs  to  be  taken  into  consideration.  Sample  cal¬ 
culations  for  circular  spreading  of  three  spill  volumes,  1000, 

10,000  and  100,000  gallons,  are  shown  in  Figures  6,  7  and  8, 
respectively.  Each  figure  shows  a  series  of  the  results  calculated 
at  a  selected  time  interval.  The  capability  of  correct  simulation 
of  the  slick  size  during  the  process  of  spreading  is  clearly  demon¬ 
strated. 
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3. 


Model  Operating  Requirements 


The  complete  computer  program  is  listed  in  Appendix  B. 

In  the  following,  the  input  information  required  by  the  program 
is  outlined. 

The  program  is  designed  for  predicting  oil  slick  movement 
in  a  harbor  or  nearshore  area.  The  primary  input,  therefore,  is 
the  information  of  the  land  boundaries.  Land  and  water  are  identified 
by  different  code  numbers.  Wherever  the  region  is  identified  as 
water,  a  water  depth  must  also  be  specified.  For  a  given  bay  area 
the  above  mentioned  data  are  constant  and  seldom  changed.  It  is 
therefore  convenient  that  these  data  are  stored  on  a  magnetic  tape 
file  in  a  prearranged  order. 

The  tidal  current  data  should  be  input  as  a  function  of  time 
and  space.  They  can  be  obtained  either  through  numerical  com¬ 
putation  from  a  separate  tidal  program  or  by  compiling  the  ob¬ 
served  information.  In  the  present  study,  a  separate  hydrodynamic 
program  is  used  for  current  data  generation.  This  program  takes 
the  same  boundary  specification  as  the  oil  slick  movement  model. 

It  requires  the  temporal  information  of  the  tidal  stage  as  input. 

The  current  magnitude  and  direction  are  computed  for  every  half 
hour  interval  and  given  as  a  function  of  location.  These  data  are 
stored  on  a  magnetic  tape  in  a  prescribed  order  and  chronological 
sequence. 

Other  information,  such  as  the  site  of  spill,  the  quantity 
of  spill,  the  oil  properties,  the  wind  speed  and  direction,  etc., 
varies  from  case  to  case.  They  must  be  updated  each  time  the 
program  is  operated. 


V. 


AN  EXAMPLE  PROBLEM 


This  section  presents  a  complete  worked  example  of  computation 
of  oil  slick  transport  using  Long  Beach  Harbor  as  a  sample  location. 

Figure  9  shows  a  map  of  the  Sail  Pedro  Bay  area.  Figure  10  shows  the 
grid-point  mesh  established  for  the  bay  and  the  water  depth  at  each  grid 
point.  The  selected  grid  spacing  is  1000  ft  and  the  asterisks  stand  for 
land  or  ocean  boundaries.  A  typical  diurnal  tide  used  for  this  study  is 
shown  in  Figure  11;  the  elevations  are  with  reference  to  a  datum  of 
mean  lower  low  water. 

Numerical  experiments  were  conducted  in  the  system  to  simulate 
spreading  and  transport  for  a  sudden  release  of  100,000  gallons  of  oil 
at  a  location  corresponding  to  grid  point  N  =  12  and  M  =  16.  Three 
experiments  were  conducted;  in  one,  only  the  tidal  current  was  responsible 
for  the  oil  transport,  while  in  the  other  two  cases  additional  effects  due 
to  wind  blowing  were  imposed.  In  each  experiment,  the  oil  release  was 
assumed  commencing  at  the  time  corresponding  hour  0  in  the  tidal 
cycle  shown  in  Figure  11.  As  discussed  previously,  a  25  hour  cycle 
tidal  current  was  generated  external  to  the  system,  but  stored  on  a 
magnetic  tape  in  sequence  of  every  half  hour  intervals,  i  Igure  12  shows 
a  series  of  plots  of  the  current  velocity  vectors  in  the  bay  calculated  at  several 
selected  reference  hours.  Similarly,  Figure  13  shows  the  plots  at  the 
same  reference  times  when  a  wind  of  5  knots  blowing  from  the  northwest 
at  284.5°  was  imposed. 

Each  experiment  was  run  for  a  period  of  36  hours  with  a  time  step 
of  6  hours.  Figure  14  shows  the  spill  site  (12,16).  Figure  15  shows  a 
typical  output  corresponding  to  18,  24,  30  and  36  hours  after  the  100,000 
gallon  oil  is  released  at  the  location  (12,16)  when  wind  is  not  taken  into 
account.  Figures  16  and  17  show  the  results  for  which  various  wind  blows 

are  included.  The  numbers  at  each  grid  point  indicate  the  thickness  of 

-2  6  2 
the  oil  film  in  10"  mm,  averaged  over  an  area  of  10  ft  around  that 

point.  The  entire  grid  space  under  these  numerals  indicate  the  area 
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covered  by  oil  on  the  water  surface  at  a  particular  lime.  The  slick 
movement  seems  extremely  sensitive  to  wind  in  the  particular  case  of 
San  Pedro  Bay;  this  is  mainly  clue  to  the  fact  that  the  tidal  currents  in 
this  area  are  generally  small. 


VI. 


FIELD  TESTS  AND  MODEL  VERIFICATION 


In  order  to  verify  the  numerical  model  there  were,  in  total,  five 
field  tests  and  demonstrations  conducted  at  San  Pedro  Bay.  The  first  two 
tests  were  made  with  3'  x  3'  cardboards  to  simulate  the  oil  transport.  In 
tne  last  three  tests,  soy  bean  oil  was  dumped  at  various  locations.  These 
tests  covered  a  v/ide  range  of  wind  speed  and  direction  and  they  were  con¬ 
ducted  at  various  stages  with  regard  to  tidal  circulations. 

1.  TEST  1,  October  5,  1973 

On  October  5,  1973  the  wind  was  principally  from  the  southern 
direction  over  the  area  of  San  Pedro  Bay.  Three  3'  x  3',  1/1  6"  thick 
posterboard  markers  were  launched  at  noon  time  in  the  southern  part 
of  the  bay  inside  the  middle  breakwater.  The  cardboards  were  sealed 
with  paraffin  on  their  edges  in  order  to  prevent  them  soaking  between 
the  layers,  and  were  painted  with  florescent  paint  of  a  bright  red-orange 
color.  Each  board  was  identified  by  a  number  painted  in  black.  The 
three  cardboards  were  launched  approximately  1000  ft.  apart  from  each 
other  from  a  25-ft  Bertram  class  motor  boat,  which  has  a  fiberglass 
hull  and  wooden  deck  finishing.  The  cardboard  markers  behaved 
extremely  well  on  the  water  surface;  they  were  flexible  and  moved 
harmonically  with  the  surface  waves.  The  positions  of  the  markers 
were  fixed  by  a  Cubic  Autotape  DM-40  electronic  positioning  system. 

The  system  includes  a  two-range  interrogator  on  board  the  boat  and 
two  responders  at  fixed  locations  on  shore.  The  Autotape  empbys 
microwave  to  measure  the  ranges  between  the  moving  boat  and  the 
two  fixed  sites  where  the  two  responders  were  located.  The  inter¬ 
rogator  visually  displays  these  ranges  once  per  second  in  5  metric 
units.  These  ranges  were  also  simultaneously  recorded  on  a  paper 
printer. 

In  the  October  5th  test,  one  responder  was  put  at  the  east  end 
of  Terminal  Island  and  the  second  one  was  at  Pier  F  of  the  Long  Beach 
Harbor.  The  exact  locations  are  shown  as  A  and  B  in  Figure  18. 
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During  the  test,  the  locations  of  the:  markers  were  traced  and  the 
distances  of  each  marker  from  the  two  responders  were  recorded 
as  a  function  of  time  (Table  1  f.  The  paths  of  the  markers  were  then 
determined  by  fixing  their  positions  using  two  range  lines  and  the 
results,  a  total  of  about  a  four-hour  tracing,  are  shown  here  in 
Figure  18.  Because  Marker  2  was  generally  following  the  same  path 
as  Marker  3,  Marker  2  was  relocated  to  a  new  location  at  2:15  PM. 

Wind  information  was  obtained  by  means  of  an  anemometer 
and  a  magnetic  compass  on  shore  where  the  responder  B  was  located. 
Wind  speed  and  direction  were  recorded  every  15  minutes;  the  records 
are  shown  in  Figure  19-  The  tidal  stage  variation  during  the  whole  day 
may  be  found  in  Tabic  II.  Taking  wind  speed  as  a  constant  of  8  knots 
but  updating  the  wind  direction  every  half  hour,  the  routes  of  the  three 
markers  were  calculated  and  compared  with  the  field  traces  as  shown 
in  Figure  20.  The  computation  used  a  500  ft  grid  space  and  a  one  hour 
time  step.  The  total  simulation  time  was  four  hours  for  Markers  1 
and  3,  and  two  hours  for  Marker  2  beginning  at  the  relocated  position. 
The  computation  considered  zero  angle  deflection  with  regard  to  the 
wind  drift  current.  The  agreement  of  the  calculated  results  with  the 
field  data  appears  to  suggest  that  the  markers  were  essentially  moving 
along  the  direction  of  the  wind. 

2.  TEST  2,  October  10,  1973 

This  was  the  second  test  to  use  cardboards  to.  simulate  oil 
movements  in  order  to  check  the  credibility  of  the  numerical  model. 

The  set-up  and  the  instruments  were  the  same  as  those  for  the  first 
test  conducted  on  October  5th,  except  that  the  wind  information  was 
read  on  the  boat  this  time.  Wind  speed  and  direction  were  read  by 
means  of  an  anemometer  and  a  magnetic  compass.  The  compass  was 
checked  against  and  in  agreement  with  the  magnetic  compass  on  board 
the  boat.  The  prevailing  wind  in  the  eariy  part  of  the  test  was  south- 
southwest;  it  turned  to  a  southwester  in  the  later  part  of  the  test.  Wind 
speed  and  direction  were  recorded  every  15  minutes  while  the  boat  made 
a  full  stop  and  the  records  are  shown  in  Figure  22. 

^Tables  are  given  at  end  of  report 
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Foul'  cardboard  markers  were  launched  inside  the  middle 
breakwater  at  approximately  noon  time,  The  four  markers  were 
placed  at  approximately  1000  ft  apart  to  simulate  a  long  slick  of 
about  3000  ft  in  length.  The  launching  points  and  the  results  of  a 
continuous  trace  of  four  hours  are  shown  in  Figure  21.  The  ranges 
of  these  markers  from  the  two  responders  A  and  B  are  given  in 
Table  III  as  a  function  of  time.  Taking  Marker  1  as  the  head  and 
Marker  4  as  the  tail  of  a  slick,  the  numerical  model  computed  a 
four  hour  movement  of  the  slick.  The  results  arc  compared  with 
the  field  data  and  show  fairly  good  agreement  (Figure  23).  The 
computation  again  used  a  500  ft  grid  space  and  a  one  hour  time 
step.  The  tidal  variation  for  the  day  may  be  found  in  Table  II.  1 
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3.  TEST  3,  February  6,  1974 

This  test  as  well  as  the  two  tests  following  differed  from 
the  previous  two  tests  in  two  aspects.  One,  real  oil  (soybean  oil) 
instead  of  cardboard  markers  was  used  in  the  field.  Two,  the 
previous  two  tests  were  essentially  to  show  the  credibility  of  the 
model  for  analysis  using  field  information;  this  test  and  the  two 
tests  following,  however,  were  to  show  the  capability  of  the  model 
for  forecasting  with  the  initial  information  and  the  capability  of 
updating  when  better  information  becomes  available. 

There  was  a  northeast  wind  over  the  San  Pedro  Bay  area 
during  most  of  the  day  on  February  6th.  The  wind  was  very  strong 
in  the  morning;  the  highest  reading  was  up  to  30  knots.  It  became 
relatively  calm  in  the  afternoon  and  the  wind  speed  was  10  knots 
at  the  time  when  oil  was  dumped. 

Instead  of  using  the  Cubic  Autotape,  a  Motorola  Mini-Ranger 
System  was  used  for  oil  slick  positioning.  Four  responders  (or 
transponders)  were  stationed  onshore;  they  were  identified  as  A, 

B,  C,  and  D  as  shown  in  Figure  24.  The  oil  dumping  and  position 
fixing  group  was  on  a. Coast  Guard  40-ft  steel-hull  cutter.  The 
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forecasting  group  was  stationed  onshore  and  equipped  with  a  portable 
teleprinter  linked  to  the  Control  Data  Kronos  Time-Sharing  System 
by  means  of  a  regular  telephone  line. 

Five  gallons  of  soybean  oilwere  dumped  at  1:33  PM.  The 
dumping  site  and  the  initial  wind  condition  v/cre  transmitted  from 
boat  to  shore  through  walkie-talkie  type  radio  communication.  The 
shore  group  immediately  fed  these  data  into  the  computer  and 
obtained  a  three-hour  prediction.  The  boat  group  continued  tracing 
the  slick  until  the  wind  stopped  blowing  at  2:29  "PM.  The  recorded 
ranges  from  the  transponders  are  given  in  Table  IV  and  the  field 
trace  is  shown  in  Figure  24.  The  predicted  results  and  the  field 
trace  are  shown  together  in  Figure  25.  The  computation  again 
adopted  a  500  ft  grid  space  and  a  one-hour  time  step;  no  angular 
deflection  with  regard  to  the  wind  drift  was  considered.  The  agree¬ 
ment  was  excellent;  the  only  drawback  was  that  the  field  trace  was 
too  short  to  compare. 

4.  TEST  4,  February  11,  1974 

This  was  the  second  field  experiment  using  soybean  oil  to 
test  the  prediction  capability  of  the  numerical  model.  The  pre¬ 
vailing  wind  was  south- southwest.  Soybean  oil  was  dumped  at 
12:15  PM  inside  the  west  end  of  the  middle  breakwater  close  to 
the  Los  Angeles  Harbor  Main  Channel.  Because  of  its  closeness 
to  the  main  channel  traffic,  the  slick  was  unfortunately  disturbed 
by  a  passing  ship,  super  tanker  ESSO  BERLIN,  and  no  accurate 
trace  could  possibly  be  made.  Despite  this  unsuccessful  start, 
a  second  dumping  was  made  at  a  new  location  at  2:13  PM.  Based 
upon  the  information  obtained  at  this  location,  the  shore  group 
predicted  the  slick  movement  for  two  hours.  The  later  field  traces 
are  shown  in  Figure  26;  the  Mini-Ranger  readings  and  the  wind 
information  are  given  in  Table  V.  The  predicted  results  together 
with  the  field  trace  are  shown  in  Figure  27. 


18 


The  tidal  stage  variation  on  February  11th  is  given  in 
Table  II.  It  is  clear  that  at  2:13  IJM  the  water  in  the  bay  was 
close  to  slack,  and  the  dominant  driving  force  was  due  to  wind 
drift.  The  predicted  results  were  obtained  based  upon  the  infor¬ 
mation  received  at  2:13  PM;  the  wind  speed  was  9  knots  and  the 
direction  was  215°  (magnetic).  As  shown  in  Figure  26,  however, 
the  field  traces  were  approximately  along  the  direction  of 
230°  -  235°  (magnetic).  The  discrepancy  was  approximately 
15-20  degrees . 

The  discrepancy  was  considered  possibly  due  to  two 
uncertainties;  (1)  the  accuracy  of  the  compass,  and  (2)  the  wind 
drift  deflection.  As  with  many  Coast  Guard  small  boats,  the  one  used 
on  February  11th  was  not  equipped  with  a  compass.  A  hand  carried 
magnetic  compass  was  used  for  wind  direction  readings;  however, 
no  error  deviation  of  the  compass  was  checked.  As  to  the  wind  drift, 
although  the  October  1973  tests  seemingly  showed  no  deflection 
effect  on  the  cardboards,  it  was  suspected  that  an  oil  slick  may 
behave  the  same  as  cardboards. 

In  order  to  clear  doubts  and  to  answer  the  above  questions, 
one  more  test  was  scheduled  on  the  following  day. 

5.  TEST  5,  February  12,  1974 

The  wind  condition  was  very  similar  to  the  previous  day. 

The  wind  direction  was  read  from  the  magnetic  compass  in  the 
same  way  as  done  previously.  In  order  to  determine  the  compass 
deviation  error  on  the  boat,  wind  direction  relative  to  a  bearing 
of  a  fixed  known  physical  object,  the  Los  Angeles  Lighthouse 
(where  the  transponder  C  was  located),  was  also  recorded;  from 
this  information  the  correct  wind  direction  in  a  magnetic  reading 
could  be  computed.  Wind  directions  obtained  by  this  method  and 
wind  directions  directly  read  from  the  magnetic  compass  are  both 
given  in  Table  VI.  Comparing  the  data  obtained  by  these  two  methods, 
one  may  find  that  on  the  boat  the  magnetic  compass  had  a  westerly 
deviation  of  approximately  10  degrees. 


19 


Five  gallons  of  soybean  oil  were  dumped  at  3:27  PM.  The 
initial  wind  speed  was  9*5  knots  and  the  correct  wind  direction  was 
220°  (magnetic).  Prediction  based  upon  this  information  again 
showed  an  approximately  10  degree  discrepancy  with  the  traced 
data  obtained  later  on.  It  appeared  that  a  10  degree  wind  drift 
deflection  had  to  be  included  in  the  computation  for  better  agree¬ 
ment.  A  complete  three-hour  trace  of  the  oil  slick  was  made; 
the  results  are  plotted  in  Figure  28.  The  predicted  oil  movement 
together  with  the  field  traces  are  shown  in  Figure  29* 

In  summary,  the  February  12th  test  led  one  to  believe  that: 

(1)  The  magnetic  compass  reading:  recorded  in  the 
'  last  three  tests  had  a  westerly  deviation  of 

10  degrees  appi oximately; 

(2)  There  was  approximately  a  10  degree  clockwise 
deflection  on  wind  induced  drift. 

Including  the  above  considerations,  a  new  prediction  of  the 
February  11th  test  was  made;  Figure  30  shows  the  predicted 
results  in  good  agreement  with  the  field  traces. 
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VII.  HARDWARE  CONSIDERATIONS 
1.  Dedicated  System 

The  desirability  of  local  operation  of  the  present  oil 
spreading  and  transport  program  prompts  consideration  of 
accommodating  the  program  for  execution  on  a  dedicated,  small 
computer.  Questions  necessarily  raised  thereby  involve  com¬ 
putational  accuracy,  execution  speed,  and  core  size  limitations. 
While  there  is  no  foreseen  difficulty  to  overcome  these  problems 
in  warranting  consideration  of  a  small  computer  configuration, 
considerable  converting  effort  for  the  present  program  would  be 
required  to  adopt  such  a  large  program  to  a  small  machine  of 
limited  storage.  In  the  following  a  proposed  hardware  configuration 
is  first  outlined,  and  then  certain  conversion  problems  are  dis¬ 
cussed. 


ll 

The  proposed  computing  system  and  its  cost  estimate  are 

i 

described  in  the  following  table. 

in 

Item 

Model 

Description 

Cost 

l  ! 

ill 

DEC  Graphic  Unit 

GT44-AA 

$34, 500 

-  16K  Processor 

PDP- 11/40  Central  Processor 

-  Scope  Display 

17" CRT  Monitor  w /Light  Pen 

!! 

-  Disk  Unit 

1.2  Megaword  Disk  Drive 
and  Control 

IT 

-  Printer 

30  cps  DEC  Writer 

11 

p 

16K  Sense  Memory 

MM  11-U 

KT  1 1-D 

Extended  memory  and 
memory  management 

7,  205 

1! 

Tape  Unit 

TM  11-FA 

7-Channel  Tape 

11, 385 

ii 

TOTAL 

$53, 090 

II 

(l 
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In  the  above  proposed  system,  the  DEC  PDF  11/40  processor 
is  adopted  as  the  baseline  machine,  around  which  peripheral  equip¬ 
ment  consisting  of  a  tape  unit,  disk,  and  display  devices  are  structured. 
It  should  be  noted  that  the  above  proposed  system  is  selected  only  as 
a  most  promising  example  and  that  it  should  not  be  considered  the 
only  configuration  satisfying  the  requirements. 

The  PDF  11/40  is  a  very  sophisticated  small  computer. 
Arithmetic  operation  and  memory  reference  times  are  approximately 
equivalent  to  those  of  the  much  larger  CDC  6600.  A  FORTRAN 
compiler  is  available  for  this  machine  and  the  DEC  assembly 
language  is  versatile  and  straight  forward.  A  magnetic  tape  unit 
is  included  for  data  handling  and  a  disk  unit  is  for  the  extension 
of  central  memory.  An  oscilloscope  display  unit  is  included  for 
interactive  program  input/output  operations. 

Based  upon  a  preliminary  analysis,  a  32K  memory  of 
the  11/40  appears  sufficient  to  contain  the  current  form  of  the 
present  program  if  single  word  integer  operations  are  utilized 
throughout  and  the  program  instruction  set  can  be  compacted 
with  no  increase  in  instruction  words.  In  its  present  form  the 
program  uses  data  arrays  at  a  prescribed  grid  scale,  appropriate 
to  an  area  equivalent  to  Long  Beach  Harbor.  If  the  present  grid 
scale  is  to  be  maintained,  any  area  larger  than  the  harbor  at 
Long  Beach  must  require  periodic  grid  displacement  to  track  the 
slick  motion.  As  the  present  scheme  stands,  the  program  updates 
the  data  arrays  every  half  hour  of  simulation  time,  corresponding 
to  the  update  intervals  of  the  tidal  current  distribution.  Each  half 
hour  of  simulation  time,  therefore,  would  initiate  a  new  compu¬ 
tational  sequence  covering  another  one  half  hour  of  simulation 
time.  Converting  of  the  present  program  would  thus  require: 
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(1)  Compacting  the  instruction  set,  via  assembly 
language,  to  reach  the  memory  capacity  of  the 
small  computer 

(2)  Incorporating  integer  arithmetic  operations 
throughout  the  program 

(3)  Developing  the  logic  to  enable  the  data  arrays 

to  follow  slick  movement  over  an  extended  area 

/ 

(4)  Developing  the  interactive  and  input /output 
capability  on  oscilloscope  display  appropriate 
for  field  application  of  the  program. 

2.  Time  Sharing  System 

The  other  candidate  for  local  operation  of  the  numerical 
model  is  a  time  sharing  system  which  permits  one  to  use  the  large 
and  sophisticated  computer,  like  the  CDC  6600,  through  a  linkage 
of  a  regular  telephone  line  to  his  local  terminal.  This  is  probably 
the  least  expensive  way  of  operation  as  the  only  hardware  required 
is  a  portable  terminal  including  a  displaying  scope.  A  software 
development  for  interactive  capabilities  on  oscilloscope  display  is 
required  as  in  the  dedicated  system.  Since  the  system  allows  one 
to  use  a  large  computer,  no  conversion  problem  is  necessary 
inasmuch  as  the  available  core  memory  is  sufficiently  large  for  the 
problem  under  consideration. 


VIII.  DISCUSSION  AND  CONCLUSION 


The  programs  presented  in  this  report  are  an  attempt  to  simulate 
by  means  of  numerical  approximations  the  movement  of  oil  slicks  in  a 
harbor  after  an  accidental  oil  spill  occurs.  This  section  reviews  some 
important  approximations,  comments  on  their  validity  and  evaluates  the 
degree  to  which  the  program  will  accurately  simulate. 

1.  Spreading  and  Transport 

The  important  approximation  used  herein  is  the  spreading 
process  simulation.  The  numerical  procedure  is  an  attempt  to 
represent  the  effect  of  the  pressure  gradient  over  the  slick  boundary 
of  the  elementary  patches.  The  concept  assumes  that  each  patch  of 
oil  slick  will  locally  spread  and  diffuse  over  a  circular  field  whose 
radius  is  stretched  approximately  by  one  grid  space  during  one 
time  step  before  the  final  area  of  the  slick  is  reached.  At  the  end 
of  each  time  step,  a  superposition  of  the  overlapping  patches  takes 
place.  It  is  obvious  that  the  grid  size  and  the  time  step  are  necessarily 
governed  by  the  spill  volume  and  the  spreading  characteristics  in  order 
that  the  spreading  process  is  properly  simulated.  This  is  equivalent 
to  the  grid  size  and  the  time  step  being  controlled  by  the  diffusion 
coefficient  in  solving  convective-diffusion  equations  by  a  finite 
difference  technique. 

While  the  necessity  of  adjusting  the  grid  size  and  the  time  step 
is  fundamental  for  a  proper  simulation  of  spreading,  the  application 
of  the  numerical  program  is  not  limited  by  these  restrictions.  From 
an  operational  point  of  view,  the  total  spreading  process  is  generally 
too  short  to  be  important.  In  other  words,  the  impact  of  slick  transport 
is  generally  more  significant  than  the  consequence  of  spreading  in 
affecting  decision  on  directing  protective  actions  and  control.  Con¬ 
sequently,  the  grid  size  and  time  step  may  not  necessarily  tie  up  to 
the  type  and  volume  of  the  oil  spill;  they  may  be  chosen  essentially 
for  convenience.  The  program  then  adopts  a  simple  procedure  to 
provide  advance  information  on  the  slick's  area  and  shape  as  a  function 
of  time. 


The  adopted  concept  of  sup  crimp  os  hip  a  transport  motion 
on  a  spreading  motion  in  calm  water  signifies  that  the  mechanism 
for  spreading  and  dispersion  is  essentially  isotropic,  which  is 
known  as  not  supported  by  observations.  Field  observations  often 
find  oil  slicks  pronouncedly  stretched  in  the  direction  of  drift. 

This  phenomenon  may  be  taken  into  consideration  in  the  numerical 
scheme  by  including  the  directional  property  in  the  spreading  coef¬ 
ficient,  provided  that  sufficiently  reliable  data  are  available  in  this 
regard. 

Field  observations  of  drift  markers  and  small  scale  oil  slick 
movements  have  essentially  validated  the  numerical  model.  While 
the  present  method  seems  to  yield  realistic  results  for  slick  move¬ 
ment  prediction,  the  more  appropriate  and  accurate  way  to  deal  with 
the  problem  would  be  to  use  a  finite  difference  scheme  to  solve  the 
diffusion  type  equation,  assuming  a  new  generation  of  larger,  faster 
computers  is  available.  At  present,  the  current  method  has  at  least 
the  following  advantages  in  comparison  with  the  finite  difference 
scheme.  First,  the  spreading  step  is  modeled  in  a  physical  way 
which  permits  immediate  comparison  with  field  observation,  so  that 
updating  of  field  results  can  be  easily  achieved.  Secondly,  the  present 
program  is  computationally  fast  which  is  not  only  important  for  quick 
decision  on  directing  control  and  cleaning  but  substantial  costs  can  be 
saved. 


2.  Sensitivity  to  Parameter  Update 

The  present  model  considers  two  primary  driving  mechanisms 
which  are  responsible  for  the  transport  of  an  oil  slick:  they  are  winds 
and  currents.  The  currents  in  San  Pedro  Bay  are  generally  small. 

The  maximum  currents  inside  the  bay  are  generally  on  the  order  of 
0.2  ft/sec.  They  are  higher  along  the  two  harbor  channels,  but 
generally  not  higher  than  0.5  ft/sec.  On  the  other  hand,  the  winds  over 
this  area  are  rather  non-uniform  and  unsteady.  As  a  result,  a  correct 
prediction  of  slick  movements  highly  depends  on  how  good  the  available 
wind  information  is.  The  computation  is  especially  sensitive  to  the 
accuracy  of  the  information  on  wind  directions. 
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The  model  is  capable  of  being  updated  with  new  information 
every  half  hour.  The  updating  may  include  any  or  all  of  the  following: 
slick  position,  wind  speed,  wird  direction,  and  current  distributions. 
Practically  speaking,  with  a  frequency  of  half  hour  updating,  the 
maximum  projection  error  in  slick  position  would  in  no  case  be 
greater  than  one  grid  space.  Theoretically,  the  resolution  depends 
upon  both  the  grid  size  and  the  total  drift  rate.  Taking  San  Pedro  Bay 
as  an  example,  assume  that  the  grid  spacing  is  1000  ft,  the  tidal 
currents  arc  negligible  and  the  average  wind  speed  is  10  knots. 

Assume  further  that  the  speed  variation  of  wind  in  a  short  time 
scale,  for  instance  a  few  hours,  would  not  exceed  +  10%  and  its 
direction  variation  would  not  exceed  jt  10  degrees.  The  average 
drift  rate  in  this  case  is  approximately  0.5ft /sec  and  an  updating 
every  three  hours  would  keep  the  projection  error  within  one  grid 
space. 


3.  Model.  Operating  Guidelines 

The  fact  that  the  computation  of  the  model  is  highly  sensitive 
to  wind  information  has  been  discussed  in  the  foregoing  paragraph. 

In  some  harbors  tidal  circulation  may  be  of  the  same  order  of 
magnitude  as  the  wind  drift;  then,  accurate  time  and  spatial  infor¬ 
mation  on  tidal  current  variations  is  also  vital  for  an  accurate 
prediction.  The  third  factor  which  may  affect  the  prediction  is  the 
information  on  the  angular  deflection  of  the  wind  induced  drift  from 
the  over- surface  wind  direction.  In  order  to  accommodate  the  model 
to  a  harbor  for  an  accurate  simulation,  matters  with  regard  to 
accurately  providing  or  collecting  the  above  three  items  must  be 
taken. into  consideration. 
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Winds  arc  known  as  the  most  sensitive  input  in  executing  the 
present  program.  To  assure  that  accurate',  temporal  wind  information 
is  available  to  the  computer,  a  number  of  remote  sensors  should  be 
established  in  the  field.  These  sensors  may  continuously  provide  wind 
speed  and  direction  information  at  various  locations  and  an  interpolation 
of  these  data  can  be  made  available  to  the  computer  at  any  time.  The 
number  of  sensors  required  is  mainly  dependent  on  the  degree  of 
uniformity  of  wind  over  the  harbor.  A  survey  by  taking  simultaneous 
recordings  of  winds  at  several  stations  around  the  harbor  can  be 
easily  conducted  to  determine  their  distribution  pattern  from  which 
the  number  and  location  of  sensors  needed  roayr  be  determined. 

Tidal  currents  are  less  critical  than  winds  as  the  current 
distributions  roughly  have  a  fixed  pattern  which  repeats  periodically. 
Problems  which  need  to  be  taken  into  consideration  are  those 
involving  systematically  taking  temporal  and  spatial  measurements. 
Field  measurements  require  substantial  undertaking  in  terms  of 
both  labor  and  time.  Computation  using  a  hydrodynamic  model  is 
probably  a  better  approach  in  most  cases  as  it  generally  provides 
better  and  more  detailed  information  than  field  data  if  the  field 
data  are  not  thorough  and  complete. 

It  is  anticipated  that  the  Coriolis  force  may  play  an  impor¬ 
tant  role  in  altering  the  course  of  an  oil  slick  if  it  is  transported 
over  a  long  distance.  Inside  a  harbor  over  a  short  time  span,  however, 
one  should  expect  the  wind  drift  deflection  to  be  negligible.  Neverthe¬ 
less,  the  field  test  of  soybean  oil  in  San  Pedro  Bay  showed  an  approxi¬ 
mately  10  degree  clockwise  deflection  in  wind  induced  drift.  There 
are  ro  sufficient  data  available  with  regard  to  the  wind  drift  deflection 
in  harbors.  Field  surveys  are  necessary  to  obtain  this  information. 
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Other  Applications 


While  the  numerical  model  was  constructed  villi  the  pur¬ 
pose  of  predicting  and  tracking  oil  slick  movement  in  harbors, 
its  application  may  not  be  limited  to  harbors  only.  The  model 
may  equally  be  valid  for  deep  ocean  and  estuaries  as  well,  pro¬ 
vided  that  some  important  considerations  are  taken  into  consid¬ 
eration.  In  the  case  of  deep  ocean  relatively  far  away  from  the 
coastal,  shoreline,  geostrophic  currents  may  be  substituted  for 
the  tidal  currents  in  the  model  to  provide  good  prediction  as  the 
tidal  excursions  are  generally  taking  an  elliptical  path  of  small 
radius  on  the  order  of  a  few  miles  and  produce  little  net  drift 
in  a  long  run.  For  application  in  an  estuary  vhere  tide  meets 
with  water  flushed  from  channels,  the  complexities  involved  arc 
to  obtain  accurate  information  of  current  distribution  around 
channel  entrances.  Current  distribution  is  a  function  of  time 
and  space.  Direct  measurements  in  the  natural  area  to  determine 
this  information  would  require  tremendous  undertaking  and  cost, 
and  generally  are  impractical.  A  hydrodynamic-numerical 
model  requires  lfss  effort  and  is  much  less  expensive,  however, 
most  existing  models  are  not  suitable  for  estuary  problems. 
Numerical  models  applied  to  oceanographic  problems  usually 
consider  them  two  dimensional.  Around  the  channel  entrances, 
however,  parameters  such  as  the  vertical  shear  distribution 
and  fluid  stratification  as  well  as  the  bottom  and  the  lateral 
friction  stresses  become  important  and  the  third  dimension 
effect  must  not  be  ignored.  In  any  case,  once  the  accurate 
current  information  becomes  available,  either  through  numerical 
computation  or  field  measurement,  the  present  model  may  be 
easily  adapted  for  use  in  estuaries. 
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A  P  P  E  N  D  I  X  A 
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Input  Data  Segue nc c  for  P r or; ranv SLICK 


Ca  rd  No. 

Colum  ns 

Va  r  iable 

Description 

1 

(8A10) 

2-80 

ID(8) 

Problem  identification 

2 

1-5 

NMAX2 

Dimension  of  computational 
grid  —  column 

(415) 

6-10 

MM  AX2 

Dimension  of  computational 
grid  —  row 

11-15 

NIND 

Total  number  of  vertical  lines 
required  to  compose  the  field 
•  points  inside  the  harbor 

16-20 

IDX 

Grid  spacing  in  feet 

3 

1-10 

Q 

Volume  of  spill,  in  gallons 

(3F1  0.  0) 

11-20 

ROIL 

3 

Density  of  oil  in  gm/cm 

21-30 

SIG 

Spreading  coefficient  in  dynes/cm 

4 

1-5 

NOY 

Column  number  indicating  the 
spill  site 

(615) 

•  6-10 

MOX 

Row  number  indicating  the 
spill  site 

11-15 

ITIME 

Reference  time  at  which  com¬ 
putation  starts,  in  hours  after 
spill  occurs 

16-20 

IMAJ  _ 

Tidal  stage  reference 

21-30 

NHOUR 

Total  number  of  hours  to  be 
computed 

5 

1-5 

IK 

0  Current  data  read  from 

Tape  7 

(15) 

'"1  Zero  currents  assumed 

7 

6-10 

IOP 

-1  Program  computes  only  the 

slick  centroid  movement,  no 
spreading  mechanism  considered 
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APPENDIX  A  (cont.  ) 


V 


Card  No,  Columns  Variable  Dose ription 

1  Both  spreading  and  transport 
are  simulated 


11-15  IT 


16-20  IP 


21-25  NP1 

26-30  NP2 

31-35  NP3 

6  1-10  WN ' 

(6 FI 0. 4)  11-20  DPH 

21-30  THM 

31-40  THMN 

41-50  THR 

51-60  ALF 


0  Ordinary  input 

1  Updating  data  required 

0  Slick  thickness  in  10  cm 
printed  in  the  output 

1  Slick  area  identified  by 

character  "O"  in  the  output 

Number  of  the  first  column  to 
be  printed 

Number  of  the  last  column  to 
be  printed 

Maximum  number  of  columns  to 
be  printed  on  each  page  of  the 
output 

Wind  speed  in  knots 

Latitude  in  degrees  at  which 
harbor  locates 

Wind  direction  in  degrees 
referred  to  magnetic  north 

Deviation  of  magnetic  north  from 
true  north,  positive  toward  east, 
degrees 

Angle  between  the  grid  column 
and  the  true  north,  degrees 

Deflection  angle  of  wind  driven 
current,  degrees  (for  zero  de¬ 
flection  angle,  a  small  number 
like  0.0001  should  be  used;  a 
zero  input  will  force  the  program 
to  calculate  the  deflection  angle 
at  every  point  according  to  the 
Ekman's  formula) 
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IF  IT  f  0,  cards. for  updating  the  slick  covering  area  are  required 


7 

1-4 

NI)XY 

Number  of  grid  points  (o  be 
updated 

(2014) 

5-8 

NN 

Column  number  of  the  grid  point 
corresponding  to  the  area  covered 
by  oil 

9-12 

MM 

Row  number  of  the  grid  point 
corresponding  to  the  area 
covered  by  oil 

(Total  NDXY  pairs  of  NN  and  MM  are 
continuously  given  in  this  card  and  the  cards 
followed) 


V 


APPENDIX  B 


PROGRAM  LISTING 
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|M>nr.[i/  m  si  ]  r; i, s  ( l :  i'i ' !  ,  mji  pH  l ,  '|  /  it  1  ,  1  A  Pi  '/  ) 
t  m  *  *  f  <  m  y  (fino  ) ,  y  r  or-f. ) ,  r  ('ion )  , nr f?on) ,  v  < <>no ) ,  ki  hivi  in  ('>oo) , 

1  u-r  i  not ,  hn>  (  i  i  oo)  ,niV‘»not ,  vw(or«o) ,  ai  pha  (<>oo)  ,iiii('»oo), 

?  V  V  f  9  o  o  1 

rfu-"OM  /mi  vi/  inno  ,  mo, i?nn  ,  i  pv 

riRMfiM  /  A I  1/  TK,  rT,  lP,MA,jnH,  1  C.  Yf.l .  F ,  lT  Yf  I  F,  JMAJ, J  \  T  Mfr  , 

*  TimM'M't 

rnMMfiM  /bay/  wi».mf  ms  i ,  mst  ,  ov ,  nvnf  hxp,ni  t  mf  m  so) , mf  ( l  so) ,  mi  (  i  so) 

rn'-'MiV!  /'II  0/  MP1  ,Mppf  MP’*,MP«»mPS 
r.iv"U'M  /ait/  M;iY,Mny,nun,sin, vr.M 

crv-Mnn  /p mi/  rnu  ,  mp,  mi  mi  ,  TOP 

BATA  1  f  Y C I  f;  ,  1  MOYPP  1/1,1/ 

(VIA  TI)TV,  J»J,  Tvf  *Ji»My,Mi<ux/lf  1  ,  1  ,o#0/ 

PI  AH  77,  in 
7  7  RV.'MAT  CSA1  0) 

i.-r a r>  so,  m,‘ax?,mhax?,mihi.), inx 
so  i  ni.'-'A r  ( i  o t s i 

T  r  f  A'B'iy ,  t  n .  0  )  NPny-MMA  X? 
ti  (nuMx.tn.  o  i  mhi’v -mm *  y? 
pr ab  si ,  o, Prni , sts 
SI  MM.'t'ATO'fl  0.0) 

RF  A  H  so,  KlnY, May,  TTT''F,TMAJ, STMT, MHnilR 
RF'AO  SO,  !Kf  inp,  TT#TP,NPl,NP?fNP3 
nniso. 

NPS=1 

I.TYfM  F  =  HHntlP/MTMT 

MAJOR=o 

Pnr:|  . 

DROr  (RD-Rnn. )  /Kfl 

nn=n 

Vri'-sO*3790. 
nyruv 
rvn=o,s*ny 
DyR=i./nx 
Ml?  1  r  MM  A  Y  ?  ♦  ^1 
mri  =MMAy?-i  n 

r.M  i  FNvinn^  cc,  Y,y,KnNvn?T,c,v,K'MAy?,M*'AY2,  i  rf  r  j nw , nr  1  ,mri  , 

*  DF  P  *  HW ,  V  W ,  Al  PHA, llll,  VV,  N«MX,  M|1»'X) 
f  Al  l  F  X  T  T 

F  Ml) 


V 


mi  P.3  i  f  \/|  |  A  i;> 


si M» t'tM 1 1  t  m  f  fiv  i  piim  ( r ,  y  ,  x,  K'lfivrc  i ,  u,  v ,  n"  a  v ty  ypf 
i  t i-'i  r;  1  ir1,  i't’1  ,  t-M-M  ,  f  f  p,mw,  vw,  a i  piia,  m 1 1  vv,  ,  mi-mx) 

in  Ml"  NS  I  MM  on*  fli"  A  v;>,  <•><  ,\  ypi ,  iiw  ( ic-1  A  yp ,  <■:*'/•>  v  p  1  ,  VW  ( R*-'A  VP,  mm  a  yp )  , 

*  Ai  i 1 » i a f nr /. v^>,  m‘<avp) 

m  i  f  n n  i nu  r  (  "max;' ,  "»•'/.  x?) ,  x ( n '* a  x  p ,  m /.  v ? ) ,  y  ( *'max?,mhax?) 

I'ltfMfilOM  Of  »(*■'  AVP,  "Ma  V,«  )  ,  vf  ni'Ay?,  mm  A  yp) 
im'i  rqni!  i'H  f  mpmx  ,  mpisv  ) ,  vv(MH"y,M‘Mix) 

rv-'i  j  s ] n n  k ri m v r  i»t  p-wav?,  wmav?y  ,  t '? f r. T i j i 1  f t « R 1  ,  mim  ) 
rm-  '-nii  /ntw/  t n  f  p  t,  *'n,i:nii  ,  1  nx 

rnt-  "fin  /ai  i  /  tk  ,  tt,  tp,  major,  i r. ycl f , mc yci  f  ,  1  ma.m  ttmf, 

*  T  O  T  V  ,  1  ,  1  M 

rnmn1  /iuy/  mt'M), nsi , MST,nx,Dyn,  nyp, mi  TMr(iso),Mr(tso),m.  (iso) 
rnp'-n"  /oil)/  ,  »‘p?,  "r>  3,  *'(>«, mps 

hoy  ,  m(iy  ,  m-m, i  r. ,  vcm 

nr  ( s?, pot ,  vcn?,?6) 


/ A»?r  / 
ropi'iin  /on/ 


rnMf'fiAj  /rr;H/  rnn  ,  *>0,  m  i  w| ,  jnp 


iniK-n  ?9,sn,3i 
31  COM  J  MUf 

on  ?o  Nr  I,  Hl’1 
nn  ?o  nr.  j ,  kp  i 
PO  t  nr  r.  T  f >N  f  M ,  M  )  =  3 
|>|P  1  ‘<  =  MR  1  -3 
WfJl  3r"Rl  -3 
nn  19  N=0,MI>13 
nn  19  m=o,yrj3 
19  iPmnNf  N,  »n  =  1 
111  I'PrHji  a  x?~? 
nn  PI  Nsl.NTMn 

HI  IMP  f  Ni)  =f-M  1 

HFfMJr? 

?\  m  |  ( M )  sr  H  m  A  V  2  "  1 
nn  in  69 
?9  PFwTnn  7 
pp  An ( 7  > 

T  f ( i vAj.pn.ft)  nn  Tn  30 
nn  ?R  1  A  =  1  r  TMA.J 
RFAnm 

?f]  rncTTMijr 

30  moon  i 

RT  A  n  f  1  )  (  fTPFr.irNf  H,  M)  ,M=1  ,MR1  )  ,Mr  1  ,  MR  ]  ) 

RPAf'dl  f  Ml  T  AT  (  M  )  ,  MF  (  K )  ,  M|  (  M  )  ,  Mr  1  ,  N  I  MH  ) 

HI  T  Mp  f  M  T  Nm  1  )=onqO 

R F  A  f>  T  1  )  (  (f)fP(M,  M)  ,  M=1  ,HMAX?)  ,M=1  ,HMAy?) 

69  PR |  Mj  70,  TO,  TDV.Mpgy 

7n  fprmatc *i*,raio, //  *  r.Rin  spat TNr.=*n o*  ft,  ttmf  spactnr 

*  13*  HH'JPS*  /) 

PRTMT  71,  MO 

7i  f np" a t ( i h° ,  *  n=*no*  r./'LKins*) 

TFfrnM.f'F.O.  )  PRTMT  IP 
ip  F'np  ^  A  T  (  i  H  1  6  X  ,  *  P F  P  hour*  ) 

PRTM1  73,  R n T I  »  S T  H 

73  Ff'PM  A  T  ( 1  HO  ,  *  PFMRTTY  MF  nil.  =  *l6.?*  (JM/FC*//  *  SPRt  ADTNH  fnf'FF 
1  F 6 .  ?  *  OYMf  S/CM*  T 


V  .... 


36 


V 


!'  (  m  .01 . 1  i  .in-1.  ,oi)  r;n  in  on 

f'Pl*'l  os 

■  ppt"1  07,  fin  iMpro),*^  f  m)  , *q_  r  w)  ,m=  i ,  mi  Mp) 

MS  f  Pl>M  A  7  (  /  0  II  MF  "(  a  ,  /  ) 

R  7  I  ni,  m»  i  m  s ) 

f'P  1  "1  VA,  (N,M=1,»'P1  ) 

oa  iimh-at  f  /  /  a  i  I'fninM  nr  sf.fiPHdM  id  rmiro  appay*//*o  m*,ootpi 
pn  0/1  *1=)  ,  M!>  i 

!’«.»  T  ►  T  o1,,  M,  M  P|  STOM (M;  >M  ,  1 1-  1  ,  KPJ  1 
/IS  |  rM  *•  A  1  f  A  A  T  ?  ,  AY  ,  0  o  1  P) 

nn  rP‘-T|>i  jr 
9 o  roin  i  i r 

[>n  r,()  ::=},mmay? 

nn  so  *■' ”  1 »  m •' a Y-> 

HrfH,")r:Vr 

so  rrs,o)  =9900. 
nn  io  wiim=  i ,  M]  nn 

M  =  M|  T  Mf  (  d!  I,M  1 
f'rr.T::*‘f  f  Ml  1-1) 

Ml  fi  T ::  i - 1  f  Ml  IM  > 

1)0  10  ’•1  =  MF f»T  ,  Ml  ST 
10  C  f  T  ,1-1  =  0. 

CAM.  MI  MDR  C  K  r  |  M  V  F  «  T  ,  M  P  A  Y  p ,  M  M  A  Y  p ,  TWIG  TOM,  Mill  ,  M  I?  \  ,  D  R  P ,  1 1 W ,  V  W ,  ALPHA  ) 

r  pn=vrM/(  i  p.  *?.s/i )  * *3 
i'n=vr'/^7oo .  i  n .  s 
1 F  f  T 1 1  110,100,110 

no  c.Aii  npnAT(r,MMAY?,MMAX?,rpn,DY,srK) 
nn  t (i  no 

100  T  F  ( P'M .  MF  .  0  1  MflrO 

r  ninivinxi  =nw  rnx*nx) 
i  '50  nn  ipn  m  =  i ,  mm  a  xp 
nn  IPO  M=1,MMAXP 

K  flMV|  P  T  f  M ,  M  )  =  1  ()  o  0  .  *  r  (  M ,  M 1  AlP.AP.SO  +0.s 
IF  (K(1NVFHT(N,M)  ,r,r.O9Q0)  f;fl  TO  IPO 
TrfJM=  KOMVK PT  f  M,  n) 

1  F ( T  T  QR , MF . 0 )  KnWVF WT(M,M)=-99 
1?0  RflL’T IfHIF 

TAM  PUT  WTC  (MOMVrPT,  MMAXP,  MMAYP,  I'P  1  ,  »  P? ,  K'P5 ,  I  TT  MF ,  HQ,  I D I  V ,  TN ,  ]  M , 
*  NRMX,VPMX) 

PALI  FIP'C  CC,Y,X,KriMV|  PT,U,V,NMAVP,MMAXP, 

1  TUfTiTHAi,  MR  1  ,MP1  ,  H  p  P ,  1 1 W ,  V  W ,  A I  PH  A ,  CPU,  i'll ,  V  V ,  S'RMX  ,  MRMX  )  ' 

PI  T  I  IP  M 
END 


V 


sm-Pf.in  1 1  ;ir  w i  vns  (w  nMvrK'  i ,  ‘may?,  '-'max?,  i lorn  him,  nwi  ,  rri  ,  pi  p,  nw,  vw, 

*  A I  n  U  A  1 

ni  mi  *'  s  i  n;i  nr  p  f  r"  ay  ?,  'm'a>v  )  ,  m  ojm/.v?,m*'ay?)  >  vu(u***y?tvvhv?\ , 

*  A  I  IMA  f  I.1'  A  V?, 

nT*‘f'f'i,'Tfvi  khmvi  i  f  * "•'  a v r? /  mmay?)  ,  t i.ti; thm (  mr <  ,  mri  ) 
r.otii'.r  im  /aii/  n  ,  n  ,  Tp.'iMmij  r.vn.  r ,  ^r.Yn.c ,  r  wa.?,  ttt  me, 

r  omt'M  / ;) a y  /  mi  »m,  hai  ,  mai  r  ny ,  dym,  nyp,  mi  t  kt ( 1  so ) ,  RF  c  1  so) ,  mi  (1  so) 
/fit  n /  mr i , 1  r ? ,  mrx » ,  mps 

WMrft. 

on  SO  M-1  , MY  AY? 
f>n  so  •  rl  ,  RMftv? 

I  ih  (  m  ,  m  )  v.  v  '■■■'  f  ri »  m  )  5  o  , 

so  rns'Mniir 

RF  AM  0)00,  MU,  MRU,  lHM,  THMM,  1 10?,  ALK 
Also  (  n  p  m  A  T ( r  F 1 0,01 
DT  M“1  I's  *■  1 1  *  I'M 
TrfWM.i-n.ft.)  on  in  ?o7 
WFPS  =  W»|*  1 . 609 
P  T  =  ^  .  1  .1 1  S° 

TMFT  A=('nMl  1  HP )  +  PT  /  1  80. 

PH  T “PPM  a  P  T /l SO . 
nn  80  Ms  t  ,  KIM  AY? 
nn  80  Mr \, MM AY? 

nf=m+? 

8  F  n  M  f  ? 

WTM'sK'r  <JS 

TF(TPTv.r,T.  1  )  nn  Tn  7^ 

T  F  (T  Pf-fi  THM  ( Ml  ,  Mr ) .  F  n .  ? )  IM- p  (  n  ,  m  )  sw  ]  md=o  . 

7  3  IF  (TK.nT.l  )  DFRfM»M)  =  10OO. 

T  F  (  A  I  F  )  7  0  #  7  S  #  7  0 
78  ALT  A=M  T*P  T  /  1  80. 

on  rn  s 

7S  f* o>?T  j  Ml  IP 

0  P  =  p  F  P  (  m  f  M ) 

FHepi  A  SORT  (A.7/ffcA.*7.?t>l>0S*SJN(i,Hn  )  ) 

APP,r?.*PI*nP/FP 
TAKA1=STWH(  A  R  P. )  -  f  >  T  M  f  A  R  (J ) 

l r ( t a m a i .FQ.n.)  no  in  ? 

T  A  M  A  ?  ~  s  f  t  j  H ( A  R  n ) f  M N  F  APQ) 

T  F  (I  A  W  A?  .  FU  ,  0  ,  )  nil  TP  3 
Al F A=A1 AM?(TAMAl ,TAMA?) 
r.n  to  s 
?  AIF  a  =  o . 
nn  td  s 

3  A  l„  r  A  s  P  T  /  ?  . 
s  rPK'T  I  'O  IF 

A I  F1 1 1 A  ( ( ! ,  MI  r  A  l.  F  A  *  1  8  0  ,  /  P  T 

|IW(M,M)r-n  .0  M(  1  .S*PI-THFT  A-Al  FA  ) 

so  vw(‘  ,m)=  n.o  ^wTMn*rns(  l  .sapi-theta-aifa) 
tf (mps.f  M.oinn  in  ?or 
?07  PPf'.'T  ?08,WMrniH,PPM,IHRf  AIF 

?08  F  n  W  R  A  T  ( F 110,  *  WIMD  S^tFO  =  *F7.a*KMnTS  THTT  A  =  *  F  7  .  ?  *  DEH  PMI  =  * 

*  F  7.?*OFP,*/3?x,  *nRT0  =*F7.?*F)En  A|.F  =  *F  7  .  ?*nEH a  ) 
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'  n  ill  m'i 
f' n <1  (  (i"i  i  *iiir 

uppn  I’lii,  WM,nTH,nI1n 

;*  1  O  f  n|;MA  1  (1  HI  ,  *  Ml  "|l  SI’f'K  I)  ::  *  f'  /  .  ?  *Y  "M  T  S  1  Ml  1  /,=  *  I  7  .  ?  *  ITO  pHl-y. 

*  i-  7.;»»iM-T.*//y  npii  t  ruifPHn  u  and  v  tm  F'T/si  r/j o<u//i 

nn  ?o?  1 '  —  i 
nn  pop 

pop  >.  omvmm  ff.'fMlrnwn  ,f)*ino. 

f-‘Di  mi  n,  (  n<nMvi  p  r  f  n,  1*1 , i ,  »j>*axp)  ,“=■  i ,  mpavp) 

rt'i'ii  in 

pn  pp7;  "=i,"max? 
nn  pr  i  hzi.umaxp 
?<va  Knt  vt  i?  r  ckj,  =vw  (  m,  *•*  1  a  i  on . 

PP  I  Ml  f  (K rir | vi  <•'  ?  Pi,  M)  ,  M-  1  ,  HeAXP)  ,  M=|  ,  M MAY?) 

PPlMl  ,>ol 

?oi  format  ( i <ii  #  *  harium<  or ptm  ju  ft  and  itptfi  angi  f  in  ornRi-TS*//) 
r*n  pon  m=i,mmaxp 
on  ?0/j  m=i,mmax? 

P  0  /I  K  niiVf  P  T  (  m,  M )  =nr  n  (  M ,  m  1 

pp  hit  v,  ( ( komvi:  "i  , «) ,  m=  i  ,  "max?)  ,  i<~  i  , mmaxpt 
ppt"1  in 

nn  prm  *'=i ,  ■•■max? 

nn  Pfi'i  h-1,"maxp 
peiG  krihVT  PI  ffj,m  =  AI  PriAfM,") 

P P  T "  1  P,  c  C KO-TVK  P  1  (  M#  Ml  ,  G  =  1  ,  MHAXP)  »  "sr  ?  ,  wu^ypi 
o  rnPMATdH  \?]U) 
io  forma t ( i ho l 

80  rn»|TlMI|r- 
PfTUPM 
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TABUS  1 


RANGE  OF  MARKER  FROM  RESPONDERS  RECORDED 
AS  A  FUNCTION  OF  TIME 


Marker  No. 


Time 

Range 

A 

(meters) 

R  ang  c 

B 

(meters ) 

Esa 

2990 

4310 

in 

2765 

4293 

13:01 

2665 

4255 

13:17 

2560 

4240 

13:33 

2450 

4196 

13:49 

2364 

4162 

14:07 

2260 

4125 

14:55 

1930 

4023 

15:20 

1821 

3996 

15:40 

1715 

3935 

15:59 

1590 

3860 

Date: 


Marker  No. 

2 

Time 

R  ang  e 
A 

(meters) 

R  ange 

B 

(meters  ) 

12:17 

3234 

4137 

12:45 

3006 

4068 

12:56 

2920 

4015 

13:09 

2780 

3975 

13:28 

2636 

3881 

13:44 

2515 

3800 

14:03 

2360 

3705 

14:12 

2290 

3660 

*14:15 

2475 

3390 

14:43 

2270 

3215 

15:06 

2120 

3100 

15:28 

2030 

2978 

15:46 

1930 

2907 

16:03 

1880 

2801 

October  5,  1973 


Marker  No. 


Time 


12:21 

12:40 

12:58 

13:14 

13:37 

13:54 

14:47 

15:12 

15:31 

15:53 


l'<  ang  e 
B 

meters  )|i met er  si 


2948 

2762 

2600 

2435 

2227 

2093 

1660 

1535 

1400 

1320 


4011 

3975 

3880 

3820 

3709 

3600 

3398 

3292 

3213 

3110 


^Relocation  of  marker 


TABLE  II 


TIDE  TABLE 
TIME  AND  HEIGHT  OF  HIGH  AND  LOW  WATERS 
IN  SAN  PEDRO  BAY 


Dale 

Time 

Height,  ft 

October  5,  1973 

0609 

3.9 

1045 

3.  1 

1614 

4.5 

2331 

0.  8 

October  10,  1973 

0126 

0.8 

0737 

5.6 

1356 

0.4 

2003 

5.0 

February  6,  1974 

1.1 

0814 

6.8 

1506 

-1.5 

2124 

4.? 

February  11,  1974 

0632 

0.  9 

1221 

3.7 

1803 

1.2 

February  12,  1974 

0047 

5.0 

0748 

1.0 

1343 

3.0 

1845 

1.8 

TAR  1,E  Ill 
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RANGE  OF  MARKER  FROM  RESPONDERS  RECORDED 
_ A S  A  FUNCT ION  OF  T I ME_ _ 

Dale:  October  10,  1973 


Marker  No. 

1 

Marker  No. 

2 

Time 

R  a.ng  c 

R  ange 

Time 

R  ange 

R  ange 

A 

B 

A 

B 

(motor  s ) 

(meters) 

(meters) 

(met  ers) 

1 1 :40 

3494 

3305 

■sa 

3373 

3487 

12:11 

3605 

2840 

mm 

3263 

3245 

12:37 

3618 

2565 

13:27 

3316 

2580 

13:48 

4227 

2002 

14:20 

3690 

2108 

14:30 

4748 

1927 

14:40 

3906 

1945 

15:25 

5398 

2040 

15:15 

4315 

1828 

15:30 

4438 

1812 

Marker  No. 

3 

1 

Marker  No. 

4 

Time 

R  ange 

Range 

Time 

Range 

Range 

A 

B 

A 

B 

(meters) 

(meters) 

(meters) 

(meters) 

ESI 

3265 

3709 

12:00 

3087 

3915 

BH 

3234 

3450 

12:30 

'3015 

3798 

13:35 

3405 

2630 

13:41 

2960 

3260 

14:22 

3760 

2080 

14:15 

3085 

2885 

14:36 

3930 

1923 

14:53 

3400 

2380 

15:10 

3550 

2180 

15:37 

3852 

1946 

53 


TABLE  IV 


WIND  SPEED  AND  DIRECTION  AND  RANGE  OF  SLICK 
CENTROID  FROM  RESPONDERS  RECORDED  AS  A 
FUNCTION  OF  TIME 


Date:  February  6,  1974 


Time 

R  ang  e 

A 

(meters) 

R  ange 

B 

(meters) 

Range 

C 

(meters) 

R  ang  e 

D 

(meters) 

Wind 
Speed 
(kts ) 

W  ind 
Direction 
(deg)* 

13:33 

2035 

3770 

10 

60 

13:58 

2240 

3510 

10 

60 

14:29 

2369 

3041 

3380 

0 

• 

*  magnetic 


TABLE  ■  V 


WIND  SPEED  AND  DIRECTION  AND  RANGE  OF  SLICK 
CENTROID  FROM  RESPONDERS  RECORDED  AS  A 
FUNCTION  OF  TIME 


Date:  February  ]],  1974 


Time 

R  ange 

A 

(meters) 

R  ange 

B 

(meters ) 

R  ange 

C 

(meters ) 

R  ang  e 

D 

(meters) 

Wind 
Speed 
(knots ) 

W  in  cl 

Dir  cction 

(dcg)v 

12:56 

3135 

4541 

871 

2215 

9-0 

215 

13:33 

3030 

4650 

770 

- 

6.  0 

207 

14:13 

2450 

3940 

1475 

1920 

9-0 

215 

14:28 

2370 

- 

1640 

2050 

11.0 

212 

14:42 

2420 

- 

1771 

2190 

11.0 

210 

14:57 

2530 

3490 

1880 

2380 

11.0 

210 

15:13 

2675 

3392 

1984 

2560 

11.0 

210 

15:26 

- 

3256 

2130 

2710 

12.0 

220 

15:43 

2760 

3120 

2290 

2840 

12.0 

208 

15:58 

2873 

2984 

2450 

3010 

12.  0 

220 

*  magnetic 
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Fieure  1  -  Slick  Area  Increased  Due  to  Surface  Tension  Spread 


Figure  5 


Nominal  Thickness  Decrease  as  a  Function  of  Time 
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Figure  6  -  Sample  Results  of  Simulation  for  Circular  Spreading  of  Oil 
on  Water  Volume  =  1000  gallons 

As  =  250  ft 
At  =1  hour 

(Numbers  indicate  film  thickness  in  10'^cm  and  numerals 
at  the  lower  right  corner  of  each  frame  indicate  time  in 
hours  after  spill  occurs.) 
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Figure  7  -  Sample  Results  of  Simulation  for  Circular  Spreading  of 
Oil  on  Water  Volume  =  1 0,  000  gallons 

As  =  500  ft 

At  =  Z  hours 

(Numbers  indicate  film  thickness  in  lO’^cm  and  numerals 
at  the  lower  right  corner  of  each  frame  indicate  time  in 
•  hours  after  spill  occurs.) 
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Figure  8  -  Sample  Results  of  Simulation  for  Circular  Spreading  of 
Oil  on  Water  Volume  =  1 00,  000  gallons 

As  =  1000  ft 
At  =  6  hours 

(Numbers  indicate  film  thickness  in  10’^cm  and  numerals 
at  the  lower  right  corner  of  each  frame  indicate  time  in 
hours  after  spill  occurs . ) 
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Figure  8  -  (continued) 
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Figure  10  -  Grid  Mesh  Used  in  the  Example  Problem  -  Grid  Spacing  As  =  1000  ft 

(Numbers  indicate  water  depth  in  feet) 


Figure  13a  -  Current  Field  at  Hour  3.5,  Wind  Speed  -  5  knots 

Direction  =  284.  5° 


VELOCITY  SCALE 


Fieure  13c  -  Current  Field  at  Hour  16,  Wind  Speed  =  5  knots 

Direction  =  284.  5° 


V 


T 


78 


iTicure  14  -  Oil  Spill  Site  at  Hour 
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Figure  15a  -  Spreading  and  Transport  of  Oil  Slick,  No  Wind 
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oil  INSTOF  t«e  bays  jooono  gallons 


Figure  15b  -  Spreading  and  Transport  of  Oil  Slick,  No  Wind 
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Figure  15c  -  Spreading  and  Transport  of  Oil  Slick,  No  Wind 
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Figure  15d  -  Spreading:  and  Transport  of  Oil  Slick,  No  Wind 
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Figure  16a  -  Spreading  and  Transport  of  Ci!  Slick,  Vr:nd  Speed  =  2  knot 

Direction  =  225° 
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icure  1 6b  -  Spreading  and  Transport  of  Oil  Slick,  Wind  Speed  =  2  knots 

Direction  =  225° 
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Figure  16c  -  Spreading  and  Transport  cf  Oil  Slick,  Wind  Speed  =  2  knots 

Direction  =  225° 
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Figure  l6d  -  Spreading  and  Transport  of  Oil  Slick,  Wind  Speed  =  2  knots 

Direction  =  225° 
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Figure  17a  -  Spreading  and  Transport  of  Oil  Slick,  Wind  Speed  =  5  knots 

Direction  =  284.5° 
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